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The superstable Weakly Imperfect Bose Gas (WIBG) was originally derived to solve the inconsistency of 
the Bogoliubov theory of superfluidity. Its grand-canonical thermodynamics was recently solved but not at 
, point of the (first order) phase transition. This paper proposes to close this gap by using the large deviations 

formalism and in particular the analysis of the Kac distribution function. It turns out that, as a function 
of the chemical potential, the discontinuity of the Bose condensate density at the phase transition point 
disappears as a function of the particle density. Indeed, the Bose condensate continuously starts at the first 
critical particle density and progressively grows but the free-energy per particle stays constant until the 
i-^ . second critical density is reached. At higher particle densities, the Bose condensate density as well as the 

free-energy per particle both increase monotonously. 

Keywords : WIBG, Bosc-Einstcin condensation, Kac distribution, large deviations, equivalence of ensembles. 



1. Introduction 

■ The proof of large deviations for the distribution of the particle density (the Kac distribution) in the Perfect 

I and in the Mean-Field Bose gases goes back to jT| . In recent papers [U [3] , the authors addressed to the large 
deviations in the particle density in a sub-domain both for the perfect and for rarified quantum gases (Fermi or 
Bose). In the present paper we extend the study of Large Deviations (LD) principle to the superstable Weakly 
Imperfect Bose Gas (WIBG) [4], known also as the Superstable Bogoliubov model [5]. The study of this model 
C3 ■ started in [H [5] was recently completed in [71 [SI [13 [H] ■ 

Actually, this model originates from a weaker truncation than that of the Bogoliubov one in the grand- 
canonical ensemble. This new system served to solve some inconsistencies between the grand-canonical Bogoli- 
ubov theory of superfluidity and the WIBG description. This non-diagonal boson model was rigorously solved 
on the thermodynamic level for the grand-canonical ensemble in [S] [TUj . It turns out that similar to the WIBG 
it manifests a phase transition with a non-conventional Bose condensation at high densities p or high inverse 
temperatures p. Meantime, even for /3 "f +oo, i.e. at a zero-temperature, only a fraction of the full density is 
in the condensate: i.e. there is a coexistence of particles inside and outside the boson condensate. This last 
phenomenon is known as depletion of the condensate. More interesting for our analysis is a discontinuity of the 
particle density from > to p+ > p_ related to a strictly positive jump of the condensate density at the 
phase transition defined by a fixed chemical potential fic- This first-order phase transition as a function of the 
chemical potential jj, sounds unusual and seems to be not quite clear as far as it concerns its physical relevance. 

In fact, the grand-canonical thermodynamics of the superstable WIBG is unknown at the point of coexistence 
of the low- and high-density phases. This paper proposes to close this gap using large deviations techniques 
description of the density distribution. For instance to answer the question, what is the value of the Bose 
condensate density when p € [yO_,p+] ? In fact several scenarios are possible. Since this phase transition is 
characterized by the appearance of a non-conventional Bose condensation, which is due to particle interaction, 
a naive thought might be that there is no condensate at all in domain p g (p_,p+), i.e. the condensate 
density jumps from zero to a strictly positive value for p > p^. In fact this scenario is wrong. Here we show 
that this discontinuity is a subtle function of the total particle density p. Formally, at the point of the phase 
coexistence, the corresponding quantum Gibbs state of the model is no more a pure state jllj but instead, a 
convex combination of some of them. A similar observation was made for example in Section 4 of [T]. 
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Actually, we verify LD for the Bose condensate density for any given particle densities p in the grand-canonical 
ensemble, i.e. even at the point of the phase transition. A direct consequence of this study is a rigorous proof 
that the discontinuity of the Bose condensate and its depletion, visible as a function of the chemical potential 
/i, does not appear in the same grand-canonical ensemble if it is considered as a function of the total particle 
density p > 0. We show that the Bose condensate density continuously increases with p > 0. In others words, 
there is no jump and the phase transition in p > is of the second order. When the particle density p (or the 
inverse temperature /3) exceeds the first critical value p_, the Bose condensate density continuously grows but 
the free-energy per particle, i.e., the corresponding chemical potential pp, stays constant: pp = pc in domain: 
p E At higher particle densities (or inverse temperatures /3), the Bose condensate as well as the 

free-energy per particle pp > pc both increase when p > p+. 

The structure of the paper is the following. In Section[2]we briefly review the grand-canonical thermodynamics 
of the superstable WIBG for a given particle density p. Our main results are formulated in Section [31 The 
proofs are collected in Section HI For the reader convenience, we collect in Appendix (Section [5]) some technical 
results as well as a short review on the LD principles. 

To conclude, we recall that throughout this paper /3 > denotes the inverse temperature, whereas p and 
p > are respectively the chemical potential and the total particle density. Also, we reserve the notation (— 
for (finite-volume) grand-canonical Gibbs state corresponding to the Hamiltonian H\. 

2. The Superstable Weakly Imperfect Bose Gas 



2.1 The Hamiltonian [4] 

Let an homogeneous gas of n spinless bosons with mass m be enclosed in a cubic box A C M"^ of volume 
V := |A|. The one-particle energy spectrum is then Sk '■= fi^k^ /2m and, using periodic boundary conditions, 
A* := [2itIj/V^/^Y is the set of wave vectors k. The considered system is with interactions defined via a 

(real) two-body soft potential ipix) ~ ip{\\x\\) such that: 

(A) Lp{x) e {E?) (absolute integrability) . 

(B) Its (real) Fourier transformation = A|n.|| satisfies: Aq > and < Afe < Ao for k G 'R^ . 

The Superstable WIBG (also known as the AVZ Hamiltonian [8] or the Superstable Bogoliubov Hamiltonian 
[7]), was proposed for the first time in j^. It is defined by 

-H'a,Ao>o '■= Ha,o + U^'^ . (2.1) 

Here the weakly imperfect Bose gas 

^A,o 2^ \ ekakttk + y ( -r^ (ofcafc + a_ka-k) + a^a.^— + —aka-k I > (2.2) 
feeA*\{o} ^ ^ ^ 

contains the kinetic-energy tern|3 plus diagonal and non-diagonal interactions. It is solved in the canonical 
ensemble in jHHO]. The repulsive interaction ensuring the superstability of by assumptions (A)-(B) is 

the "forward scattering" interaction 

C/f := ^ E = ^ (^A - ' ^i^^ ^A := E "fe"'^ (2-3) 

fci,fc2eA* fcSA* 

defined as the particle-number operator within the grand-canonical framework. Indeed, aj! and Uk are the usual 
boson creation / annihilation operators in the one-particle statclll XA(a;)e*'^'^/v^, acting on the boson Fock space 

:=0 n^B \ with n^-^ := [L' m)„ , C, (2.4) 

n=Q 

* Recall that eq = 0. 

tHere XA i^) is the characteristic function of the box A. 
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defined as the symmetrized ??-particle Hilbert spaces, see [TTl [12] . 

Remark 2.1 Let Hqa C (A) be the one- dimensional subspace generated by tpk=o (x) = 1/VV. Then « 
J-QA ® J- A where J-qa and J-j^ are the boson Fock spaces constructed out of TioA o,nd of its orthogonal complement 
Ti-QA respectively. 



2.2 Grand-canonical thermodynamics for a fixed particle density [HI [9l 110] 

We consider liere tlie grand-canonical ensemble ^) defined by a given particle density p, or more precisely 
by the chemical potential )J.a,p-, which is a unique solution of the equation (|2.5p below. In any finite volume, 
the corresponding particle density is strictly increasing by strict convexity of the pressure. Therefore, for any 
p > 0, there exists a unique fiA,p such that 



(2.5) 



A,Ao 



where (— )^sb always represents the (finite volume) grand-canonical Gibbs states for H^^^ taken at inverse 



temperature (3 and chemical potential ^A,p. In the thermodynamic limit, ^A,p converges to G IB 
p > 0. In fact, ^p is strictly increasing except for p G where it equals fXc = ^ic- Here p+ > p_ 



for any 
> are 



two well-defined density only depending on the inverse temperature /3 > 0. Additionally, /ip := Up + AqP with 
Q!p < and d\„ap = for p e [p-,p+]. 

Moreover, there is a non-conventional Bose condensation induced by the non-diagonal interaction Uj^^ for 
high particle densities: 



=lim 

A 



V 



H 



A.Ao 



= for /9 < p_ , 
> for p > p+, 



(2.6) 



with d\gXp = for p ^ [/9_, p+J. When p i note that the Bose condensate density Xp converges to Xp^ > 0. 
In particular, since pp = pc for p G the Bose condensate density a;^ as a function of the chemical 

potential p jumps from to Xp^ at p ~ pc- An illustration of the behavior of Xp for a fixed density p (or a;^ at 
a fixed chemical potential p) is performed in Figure [TJ 





To zero-temperature 



To zero-temperature 



Figure 1: Illustration of the Bose condensate density, Xp at fixed particle density p > or x^^ at fixed chemical 
potential p £ R. The dashed line closing continuously the gap between p_ and pj^ in the illustration of Xp 
is a consequence of results of the present paper. Here each of the asymptotic straight lines are : Xp = p, or 



Xp = p/\o. They correspond to the limits : Xp 



or Xp^oo) with 100% of the Bose condensate. 



We would like to stress that coexistence of different types of condensations is a subtle matter. For example a 
slight modification of interaction, see [H [Q] [10] , excludes any coexistence of non-conventional and conventional 
Bose condensation, as it appears for high densities in the Bogoliubov WIBG [T51 114j . 
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Below we consider coexistence a high and low density phases. For intermediate total density p (f. p+] we 
find for the particle density outside the zero-mode 

where fk eu - ap + XpA^ and Ek := (/^ - 2:^A^)^. 

Observe that the grand-canonical thermodynamic behavior of the superstable WIBG is unknown for p e 
at fixed /3 > 0. When (3 +oo, i.e. at zero temperature, the critical densities p_ and p+ could 
both converge to zero depending on the interaction potential O [10] , whereas the critical chemical potential jic 
converges to a negative value. Moreover, we have a non-zero particle density outside the zero-mode for any fixed 
p > even at zero-temperature since 

lim lim {alak)irSB > and lim hm ( \ < p. (2.8) 

/3->+oo A V , ^'^0 /3-^+oo A \ V „SB 

In other words, there is a depletion of the Bose condensate even at zero temperature. 

To conclude, the grand-canonical pressure associated with H^^^ in the thermodynamic limit equals 

= supjinf (po^(/3,a,a;)-f ^^^^-^11 (2.9) 
ri->0 L"<o L ^Ao J J 

= M \p^iP,a,x,)+^-i^^^-^] (2.10) 

- p^(/3,ap,xp) + ^p^ (2.11) 

for any p > 0. In particular, Xp and Up are solutions of the variational problems (|2.9p and (|2.10[) respectively. 
For any x > and a < 0, 

(/3, a, x) :=lim ^ InTr^; |g-/3(i/- „(.,«)-«.) | ^2.12) 
is here the (infinite volume) pressure of the so-called Bogoliubov approximatiorQ 

H^Q{x,a):= ^ [{eu - a) alau + {alak + a_f.a^k + ala_t. + aka-k)\ (2.13) 

of {H^Q—a{NA—a'^aQ)}. Observe that H^q {x, a) is defined on the boson Fock space J^'j^ for non-zero momentum 
bosons, cf. Remark 12.11 and Section [3] for a rigorous definition of the so-called Bogoliubov approximation. 
Finally, note that the Hamiltonian ^ {x, a) represents, via a unitary transformation, a perfect Bose gas of 
quasi-particles with one-particle spectrum Ek for k E A*\ {0}, see for instance [5]. 

3. Large deviations for the Bose condensate at a fixed total particle density 

To define the (finite volume) distribution Da_p of the condensate, we first recall the rigorous definition of 
the Bogoliubov approximation due to Ginibre [15j and based on coherent vectors. For any complex c € C, a 
coherent vector |c) is an element of the boson Fock space JFqa for zero momentum bosons (cf. Remark 12. ip . 
satisfying ao\c) = c^/V\c). In fact, if Qq is the vacuum of JF|*, then |c) := exp{— F|cp/2 + cy/Va'^}i}o for any 
c € C. The Bogoliubov approximation of a self-adjoint operator A acting on JF^ is the operator A(c) defined 
on the boson Fock space J^'j^ without the zero mode by its quadratic form 

(VilA(c)IV^) :=(c®V'i|A|c®V2), (3.1) 



■t Combined with a gauge transformation aj^ 
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for |c ® 2) ™ the form-domain of A. 

Now, for any chemical potential /i G M the (finite volume) grand-canonical pressure associated with -fff Aq 
equals 

(/?, := ^ InTr^^- {W^} , with Wa := e-'^f^A^o-^^-). (3.2) 

By using the generating family of coherent vectors |c) for c G C, we can rewrite the trace Tr above to observe 
that 



(/3,M) = ^l-i / TrniW^A W}d^c= ^Ini- J e^W^^M^^ (3.3) 



where d^c := ^dcidc2 with c := ci + ic2, Wa(c) results from the Bogoliubov approximation (|3.ip of the 
statistical operator Wa, and 

pr(/3,/.,c):=^lnTr^, {W^a(c)} (3.4) 

is the pressure defined by the partial trace. For any p > 0, the corresponding distribution Da./^ related to the 
Bose condensate number density, is now defined on the Borel subsets ^ of C by 

Ba,^ [A] := e-'^^fA"^^'^)^ j e'3^fi"('3''''=)d2c. (3.5) 

A 

Then, at fixed particle density p > 0, we express a large deviations principle (Section 15. 2p for the condensate 
distribution Da,p := '^A,^^A.p■ 



Theorem 3.1 (LD principle for the condensate distribution at a fixed density p) 

The sequence {Da.p} satisfies a large deviation principle with speed PV and rate function 



Dp (x) := irif < pg {P, a, Xp) + — > - inf <^ Pq {(3, a, x) + 



for X = |c|2 > 0, c/. fOD- iIXT^] . 



This theorem shows in particular that the probability to observe a density uq/V G ^ of condensed bosons 
enclosed in A for a fixed chemical potential G M decreases exponentially with the volume V = |A| if the 
distance between the Bose condensate density Xp (|2.6p and the set C M is strictly positive. Now, the next 
step is to evaluate the limiting probability measure, in particular at the phase transition defined for a chemical 
potential p G Recall that the Bose condensate density Xp (|2.6p converges to when p ] p- but to a 

strictly positive value Xp^ > when p I p+. 



Theorem 3.2 (The condensate distribution outside the point of the phase transition) 

The (finite volume) distribution ©A.p of the condensate converges weakly in the set of probability measures 
Ml (C) as A 1 M.^ towards the singular measures 

Dp :=lim Da,p j 5(c^ xl'^e'") AO, 



for any p G (0, p_) U (p+, -l-oo). 



For /3 -t-00, i.e. at zero-temperature, observe that p_ and p_|_ could both converge to zero, depending on 
the interaction potential. But. at finite temperature, i.e. at /3 > 0, one always has p+ > p_ and the convergence 
of Da.p is not solved for p G [p-, p+]. The corresponding result is therefore expressed in the next theorem. 
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Theorem 3.3 (The condensate distribution at the point of the phase transition) 

Let p+ > As A I M'^, the (finite volume) distribution ©A.p of the condensate converges weakly in A^i (C) 
towards a convex combination of the singular measures 



2ir 



:=lim ]D)A,p = (l-^^)5{c) + ^ j 5{c- xjf e^«) d0, 



for any p G [p^ , p+] and with Kp := (p — p^)/{p^ ~ P-)- 

Note that Kp is a strictly increasing and continuous function from [/?_, p+J to [0, 1] . This rcsuh gives a strong 
evidence that, at the phase transition, the corresponding Gibbs state is not a pure state anymore but a 
convex combination of pure states, see for example Section 4 in [T]. 

Integrating Da^p with the function ip{c) = |cp, we finally obtain the Bose condensate density (|2.6p inside the 
phase transition, i.e. for p G [p^,p^]. 

Corollary 3.4 (Derivation of the Bose condensate density for any total density) 

The Bose condensate density equals 

{0 forp<p_. 
— — —Xp^ for p e 
P+-P- 
Xp>0 foTp>p+. 

In particular, it is continuous as a function of p > and linearly increasing for p G [p_ , p+J , cf. figure QJ 

As a function of the density p > in the grand-canonical ensemble, the phase transition is of order two if 
p+ > p- whereas it is of order one as a function of the chemical potential. In particular, take p < p_, then the 
system behaves as the so-called Mean-Field Bose Gas, i.e. the model defined by the Hamiltonian 

E ^^"^^'^ + ^ - A^a) , (3.6) 

with no Bose condensations. Increase now the particle density. The free-energy per particle, i.e., the chemical 
potential pp^p < pc, normally grows until we reach p = p^. By further increasing of the density, a Bose 
condensation continuously appears to reach the value Xp^ for p = p^. Meanwhile, the corresponding chemical 
potential pp stays constant at the phase transition: pp = pc for p € [p_, p+]. Finally, at higher particle densities, 
i.e., for p > p+, the Bose condensate as well as the free-energy per particle pp > pc both increase. 

4. Proofs: Large Deviations for a generalized Kac distribution 

We are going to study the grand-canonical ensemble at a fixed total particle density p > 0. But before doing 
this, we start our analysis at a fixed chemical potential p. Then we prove the LD principle for the condensate 
plus "out of condensate" particle densities. The corresponding distribution Ka,/j is a combination of the so- 
called Kac distribution [T] for particles outside the condensate with the condensate distribution Da,;^. This is 
expressed by Theorem 14. 1[ which is therefore, a generalization of Theorem 13.11 To study the phase transition, 
we use the generalized quasi-average procedure [T] by taking a "perturbed" chemical potential 

Ac := A^c + ^ + o (^^) for7eM, (4.1) 

we analyze the thermodynamic limit of the generalized Kac distribution at this chemical potential, sec Theorem 
14:21 



As a consequence, the generalized quasi-average procedure (|4.1[) gives the finite volume behavior of the 
chemical potential pA.p solution of (j2.5|) at the phase transition, i.e. when p e [p_,p+] if p_|- > p_. Indeed, by 
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applying the distribution Ka,^ to an appropriate function, we obtain the mean particle density at a chemical 
potential flc for any 7 € M. This procedure will then imply that for p £ there is a unique and explicit 

7p such that /iA,p = fic with \jp\ = o(y), see Section |4?2] 

Meanwhile, the large deviation principle for Ka,^ given by Theorem 14. II directlv implies Theorem 13. II for any 
p > 0. Applying the result of Theorem 14.21 to the chemical potential /iA,p = Mc for 1 = Ip, we also get Theorem 
13.31 for p S [p-,p+\. If p ^ [p-,p+], the generalized quasi-average procedure is not necessary and Theorem 13.21 
is a simple consequence of Theorem 13.11 We give now the promised proofs. 

4.1 Large deviations for generalized Kac distribution 

The particle number density as a M-valued random variable, is well-defined via a well-known probability 
measure, the so-called Kac distribution [T]. We give here a generalized version of the Kac distribution associated 
with the condensate and its depiction. This distribution is defined, on the Borcl subsets ^ C C and B C M.^ by 
integration over the zero-mode coherent state: 

Ka,p [A] e-^^P-"^^^''^^ I d'c I vj, (dy) e^^(K^+I^O-/f"(/3-^.^)), (4.2) 

A B 

with 

+00 

i^A(dy) :=^,5([yF]-n)d2/. (4.3) 

n=l 

Here [.] is the integer part and 

/f«(/3,y,c) := --l,lnTr^,v,^^ [{W^ (,)}([.v],.^o)^) (4 4^ 

where Wa,o(c) results from the Bogoliubov approximation (|3.ip of the statistical operator Wa,o P-2p and 
^(n,A:#o) ^jjg restriction of any operator A acting on the boson Fock space T'g to the space {'^qa}^"'' of 
non-zero momentum bosons. Now we express our first result concerning large deviations for the generalized 
Kac distribution Ka,^. 

Theorem 4.1 (LD principle for the generaUzed Kac distribution) 

The sequence {Ka,^} satisfies a large deviation principle with speed f3V and rate function 

K^ (x, y) (/?, p) + f^ (/?, y,x) + ^ [y + xf - p [y + x) . 

Here x = |cp > 0, y > and 

fo {P,y,x) :=sup {a{y + x) - p^ {P,a,x)} 

a<Q 

is the Legendre-Fenchel transform of Pq {f3,a,x) i2.1^) . 

Proof. Let us start by some observations. The pressure Pq {P,a,x) defined in (|2.12p and used in (|2.9p can be 
explicitly computed. Indeed, 

p^{(3,a,x) = ax ^ J In (l - e-'3V(^'=-")(-''^-"+2^^'=)^ d^fc + 

+ / ^3 / {ek-a + x\k - Vi^k -a){ek-a + 2x\k)} d^k, (4.5) 
2 (27r) J ^ -* 

R3 

for any a < 0. Since p^ (/3, a, x) is a convex function of a < 0, it is also the Legendre-Fenchel transform of 
fo y, , i.e. 

Po ifi, a, x) =sup {a{y + x) - f^ {(3, y, x)} for any a < 0. (4.6) 

y>o 
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Combined with (|2.9[) this last incquahty imphcs that 



(/3, m) -sup ( inf (sup (a(y + x) - {P, y, x) + ]]] ■ (4-7) 

We would like to bring the infimum over a < inside the two other suprema. In general, a supremum and 
an infimum do not commute. In this peculiar case, this is however the case. Indeed, for any fixed a; > the 
function ^ 

^ [y, a) a[y + x) (/3, y, x) + (4.8) 

is a strictly concave function of y > and a strictly convex function of a < 0. Then, we obtain the uniqueness 
of the stationary point (j/, a) solution of 

dy^ {y, a) = and c)^* {y,a) = y + x+ ——!- = 0. (4.9) 

In particular, we can commute the infimum over a < and the supremum over y > in (|4.7p to obtain 

/^(/3,m)= sup \fi{y + x)~fo''i(3,y,x)^^iy + xf]. (4.10) 

This result is coherent with the rate function ix,y). By explicit computations, observe also that there are 
M, B > such that any solution (x^, y^) of the variational problem ()4.10|) verifies < M and j/^ < M whereas 
for any x > M and y > M we have 

M (y + a:) - (/3, y,x)~^{y + xf <-B{y + x). (4.11) 



Now we are in position to analyze the LD principle for distribution Ka,;^ (Section 15. 2p . 

From (j4.5p the rate function (.t, y) is not identical oo and has compact level sets, i.e. for each m < oo, the 
subset {{x, y) : (a;, y) < m} is compact. 

Let a closed set C Co x Ci C C x R_(_. Remark that M can be taken arbitrary large (and B being the 
same). Then, without lost of generality, we can assume that any c G Co and y € Ci satisfy |cp < M and y < M 
respectively. By (|4.1ip . we also obtain 

Co Ci 

_^l.^^^PV{priP,f.)+2BM} J J (dy) ^^-pVB(\cr-+y)_ ^^^^^^ 

C M+ 

For large enough M, one has 

2BM+ sup l^fiiy + x) ~ fi {(3,y,x) - ^ {y + xf"^ > 0. (4.13) 

Consequently, the inequality (|4.12p combined Lemma [5.21 implies that 

limsup -^IuKa,^ [C] < - inf {\c\^y) . (4.14) 

In other words, the large deviations upper bound (|5.22[) for Ka^^ with speed /3V and rate function is verified. 
It remains to analyze the corresponding large deviations lower bound (|5.23p . 

Let G be an arbitrary open subset of C x R_|_. Note that 

Ka.^ [G] > Ka,^ [{(c, y)}] - e-/5^Pr(/3,M)e/3v{M([yV]+|cP)-/f «(0,,,c)}^ 
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with (c, y) G G- From Lemma |5.2[ it yields that 



^ InKA,^ [G] > -K^ {\c\\y) . (4.16) 

Since the last inequality holds for each point of CJ, it means that 

hm^inf ^ InKA,^ [G] > - inf {\c\^ y) , (4.17) 

i.e. the corresponding large deviation lower bound (|5.23p for Ka ^ holds with speed f3V and rate function K^. 
□ 

For II ^ fic 'Vfe already know that the variational problem (j4.10[) has a unique solution (a;^,?/^). Therefore, 
as a direct consequence of the fact that the sequence {Ka,^} satisfies a large deviations principle with rate 
function having a unique minimum in at {x^,y^) for any n ^ /itc, the distribution Ka,^ converges weakly 
on the set of probability measures A^i (C x M_,_) as A "f M'^ towards the singular measure 



=lim Ka,^ = ^ J ^i''- ^ - d'^' ^' + (4.18) 



Now, the next step is to evaluate the limiting probability measure at the phase transition defined for a chemical 
potential fi = fXc- Indeed, if p+ > p-, the solution (x^,?/^) jumps when p, cross the critical chemical potential 
He from (0, p_) to (xp^, yp+) with Xp^ > and yp_^ := p+ - Xp^ > p_. 

Theorem 4.2 (The generalized Kac distribution at the phase transition) 

If p+ > then the distribution ^A,fi^ converges weakly in Mi {C x R_|_) as A "f R'^ towards 

hm IKa,p. = ^^S (c) <5 (y - P-) + J S[c- x],{^e^') S {y - y^J dO, (4.19) 



with := (1 + ef^P+-P-'^)-^ G (0, 1) and fi^ defined by (gjj) for any 7 € R. 

Proof. We have already mentioned that the rate function K^^ has two distinct minima in R^ at (0, p_) and 
{xpj^,ypj^)- To get around this complication, take e G (0, Xp^) D (0,?/p^ — p_) and define 

A- := {ceC: |c|2 G (0,Xp^ -e]} x (p_,yp^ -e] (4.20) 

and 

^+ := {cG C : |cp G (xp_, -e,+oo)} x (yp^ - e, +00) . (4.21) 

Now, let K, - and K„+ be defined as the two restrictions of Ku to A- and A+ respectively and remark that K, - 
and have both a unique minimizer in R'^, respectively (0,p_) and {xp^,yp^). Define the corresponding 
probability measures 



which satisfy a large deviations principle respectively with rate functions K^- and (where x = |cp). Take 
any positive and continuous function Lp{c, y) of (c, y) G C x R+ and observe that 

/d^c / z.A(d2/)v:>(c,y)e'^^('^=(^+l=l')--^'"(^''y^=') 

/ if (c, y) Ka u A (d^c) = -, ^ $X + (4-23) 

i/^ '^^"^ ^ /d^c / z.A(dy)e''^('^=(^+l^l='-^A"('''^^=') ^ ^ 
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with 

■ " / e{7+o(i)}(j/+|cP)L- (d2cdy) + Oa / e{'r+o(i)}(y+kP)L+ (d2cdy) ' ^^'^^^ 

/ (^(c,2;)e^^+°(i)>(^+l^l')L+ (d^cdy) 

■ " / e{'/+°(i)}(y+l=l')LX(d2cdy)+ / e{7+o(i)}(?/+|cP)L+ (d2cdy) ' ^^'^^^ 

J g/3V{Mc(«+|cn-/r(/3,.,c)}^^ (d2/)d2c 

®^ J e/3^{^.to+kP)-/f-(3,y,c)}^^ (dy)d2c' ^^'^^^ 

By Lemma [5.21 the fmiction 

/^c(2/+|cn-/r(/3,2/,c) (4.27) 

converges in the thermodynamic hmit to 

f^ciy + Icp) - /o^ (/?, y, |cp) - ^ (y + \c\'f , (4.28) 



and 



which has suprema at (0, p-) and {e^^ Xp^,yp^) for any 6 € [0, 27r]. Consequently, the coefficient 0a (j4.26p 
converges to 1 in the thermodynamic hmit. Since p+ — yp_^ + Xp^, it is then straightforward to see that 



2ir 



hm <&X = ^7^(0, P-) and lim $+ = 1 v'(.Ty/ y,Jd0. (4.29) 



Let us apply these limits to the function ip{c,y) = e~*^'^+''^ with t > 0. Then, by bijectivity of the Laplace 
transform, it follows that IKa./i^ converges weakly on A^i (C x R+) as A "f R'^ to (|4.19p . □ 

Notice that the function : R ^ (0, 1) defined in Theorem 14.21 is strictly decreasing and in fact bijective. 
Therefore, by applying Ka^/j^ to (p(c, y) ~ |cp + y, we have shown that the particle density can converge to any 
fixed density in the open set (p_, p+) : 



Aa\ 



hm^-^^^^ =CyP- + {l-^^)p+- (4.30) 



-"A, An 



Note that all these results are coherent since we have 

p_ = lim {^^p- + (1 - £,jP+)} and p+ = lim {^.yP- + (1 - £.jP+)} ■ (4.31) 

7 — > — oo 7 — »+oo 

In particular, if 7 = 7a = o{±V) in (|4.ip diverges to ±00, then we would obtain one of the previous limit, 
depending if 7a i —00 or 7a T +00. 



4.2 Application of the generalized Kac distribution for a fixed particle density 

Let us consider now the total particle density as a parameter that defines the grand-canonical ensemble. 
Theorems 13.11 and 13.31 are direct consequences respectively of Theorem 14. 1 1 and (|4.18p for the chemical potential 
Hp defined as the thermodynamic limit of fi^.p (|2.5p . The only remaining question is to study the case of fixed 
particle densities at the the point of phase transition, i.e. in domain: p G (/9_,p+) for p+ > p_. From (|4.30p . 
we obtain that 

lim/^\ =pe(p_,p+), (4.32) 
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for a chemical potential 



cf. (HH). Therefore, 



^ + ^ + (l^) ^itl^ > •= ) , (4.33) 

pv \f3vj P+-P- \p+-p; 



MA.-Mc + |. + o(-]. (4.34) 

In particular, from (|4.2[) with j = jp wc get Theorem 13. 31 for p G {p-,p+). Recall also (|4.3ip . In other words, if 
p = then 7p < (|7p| = o(y)) would diverges to — oo, whereas if p = then jp = o (V) +oo. It follows 
that Theorem 13.31 is proven for any p G 

5. Appendix 

In this appendix, we first give supplementary results needed in the previous section. Next, for the convenience 
of our reader, we shortly repeat the notion of large deviations principles. 

5.1 Some technical statements and proofs 

The thermodynamic limit of pf^(/3,p,c) (|3.4p is first analyzed in order to obtain next the one of the free- 
energy density fl^{(3, y, c) (|4.4p . which is given in Lemma [5T2l 



Lemma 5.1 (The pressure p^^(/3, /i;c) in the thermodynamic limit) 

For any c G C, /i G K and (3 > 0, the pressure p^^{P, p, c) converges towards 

/^(/3, p, c) :=lim pi^(/3, p, c) = inf \p^{(3, a, x) + ' 



A a<0 ZAq 

Here x = |cp > and recall that p^ {(3, a, x) is defined in 112.12]) . cf. also 



Proof. The proof is obtained by a comparison between suitable lower and upper bounds for pf^(/?, /i, c). We 
start by the lower bound. By taking any othonormal basis {(V'nl}^! of T'^, 

oo 

Tr^, {Wa (c)} =J2 (c® <l e-'^^^-'^o-^^A) |c® , (5.1) 

n=l 

and so, by the Peierls-Bogoliubov inequality we get 

Tr^; {Wa (c)} > sup jfj ^-/'(-^'A^l^r.o-A'^Alc^V.;) | ^ Tr^, |g-/3ffIio(c,M)| ^ (5.2) 

see e.g. [IllIlT], where H^^^{c, p) results from the Bogoliubov approximation (j3.H) of — pNa}. From 

[5] we already know that 



lim ^InTr^, [e~^"'^^o(-^)} =inf |po^(/3,a, |cp) + 



(m - a)^ 
2Ao 

Consequently, the inequality (j5.2p implies in the thermodynamic limit the lower bound 



(5.3) 



/^(/3, /i, c) > inf \ p^iP, a, |c|2) + ^^^^^ ^ , (5.4) 



for any c G C, p G M and /3 > 0. 
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To obtain an upper bound on p^^ {(3,ii,c), wc follow the idea of [TBI , and use the coherent state representation 
of {^A,Ao - /^^a} given by 

H^H,, -fiNA^J d'c {h^X: m) |c) } ' (5-5) 

c 

where the Hamiltonian H^^^ (c, /x) is defined on J^'j^ by 

^A,l(c,/i):=i?Al(c,M) + A, (5-6) 

with 

V V 



A:=/.-2Ao|cp + ^-l ^ (Ao + Afc)4a;.-. (5.7) 



fceA*\{o} 



Actually, -ff^ Ao ("-'M) derived by replacing the operators a|5ao, aoao, ajj'^O' '^o'^o'^oflo in {-^a Ao ^ M-^a} 
respectively by |T^cp — 1, Vc^, Vc^ and y^lcl"* — 4y|cp + 2. Let {(?/)^ (c) be an othonormal basis of 

eigenvectors of Ao ('^' ■ Since for any z, c G C 



it follows that 



-pfd^zH^I (z,p)|z)(z| 

Tr^;{WA(c)} = ^(c®VUc)|e - Ic^VUc)) 

m=l 

X n « ^A,Ao (^J' M) K (c)) \ , (5.9) 

with the two real- valued functions Rm and I„i of (^i, • • • , z^) S C™ defined by 

2 ™ 2 
R™ (zi, • • • ,z™) := |zi - c| + kj-i ~ ^jT + km ~ c| , 

Im (zi, • • • , z„) := i (zic - CZi) + i X) (^j-^j-i " ^j-i^i) + * (c^m - ZmC) . 



(5.10) 



Since Im (c, • • • , c) = and 



inf Rm (zi, • • • , Zm) = R„i (c, • • • , c) = 0, (5.11) 

(Zi,--- ,Zm)eC™ 

by virtue of (|5.9p combined with large deviations arguments, one can obtain in the thermodynamic limit that 

p^^(/3,M,c) =lim ^InTr^; je'^^rAo^^^^) } . (5.12) 

Justification of the LD technique in sums (j5.9p is based on the uniform domination theorem and it follows the 
line of reasoning developed in (TQ. Meanwhile, by using the Bogoliubov convexity inequality [M] it follows that 



Tr,, {e~^^n.i^-^] < 2_i,Tr,, ,io(-) } - ^ (A)^.,^^,^^^ , (5.13) 



where 

Tr^, {-e-'^^A'of-'^)} 



(5.14) 
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In particular, since for fc S R"^, < < Aq by our assumption (B) on the interaction potential, we obtain from 
the inequality (|5.13p together with (|5.7p that 



. In Tr^, {e-'^^A-o t^-'^) } < ^ In Tr^, {e-^^A.'o } + 



/3y -"A I J - -^-^ r J V T/2 

A:GA*\{0} ' ° 

The last term can be explicitly computed. We omit the details. In fact, for any /i e M one can check that 

^ E K«fe>ifs.^(c,^)=0(l) asATK^ (5.16) 
feeA*\{o} ■ ° 

Therefore, from ()5.15p together with ()5.3p and ()5.12p one deduces that 



p^^(/3,/i,c) <ini\p^{(3,a, \c\^) + '-^^-^ \ . (5.17) 



(m^)2 

--<o [^^> v-'-'i-i / - 2Ao 
Together with the lower bound (j5.4[) . this inequality proves the lemma. □ 
Lemma 5.2 (The free-energy density f^^{P,y,c) in the thermodynamic limit) 

For any c G C, y > and /3 > 0, the thermodynamic limit f^^{(3,y,c) of the free-energy density f^^{(3,y,c) 
Ii4-4\ l equals 

/^^(/3, y, c) :=lim /f ^(/3, y, c) = f^ (/3, y,x) + ^{y + xf , 

wii/i X = |cp > 0, and f^ (/3, y, a;) defined as the Legendre-Fenchel transform of Pq (/?, a, x) i2.12^) . cf. Theorem 
1\ 



Proof. The pressure pf^^{f3, fi, c) (|3.4p can be rewritten as 



Pr(/3,M,c) = ^ln 1 e''^(''^-/A"(^^«-^))^A(dy)+M|cp, (5.18) 

M+ 

with i/A(d?/) defined in (|4.3p . It is then straightforward to check that the thermodynamic limit p^^ {f3, iJ.,c) of 
pf^(/3,^,c) (EH) equals 

/^(/3, ^, c) =sup {^y - /^"^ (/?, y, c)} + /i|cp, (5.19) 

with f^^{f3,y,c) < oo for y > 0. The derivative of the pressure p^^ {(3, iJ,,c) is continuous as a function of /x, 
cf. Lemma 15.11 and (|4.5p . Thus, by using the Tauherien theorem proven in I20|. the existence of p^^ {(3, ij,,c) 
already implies the convexity of f^^{f3,y,c) as a function of ?/ > 0. In particular, it yields that 

/^^(/3,y,c)=sup{M(j/ + |cn-/^(/3,M,c)} for y > 0. (5.20) 
By using the explicit form of p^^{P, p,, c) given by Lemma 1 5. 11 a straightforward computation then gives 

/^^ (/3, y, c) =sup {a(2/ + \cf) - pi{(3, a, \c\^)} + ^ {y + x)\ (5.21) 

a<0 ^ 

□ 

5.2 Large deviations principles 

Let X denote a topological vector space. A lower semi-continuous function I : X ^ [0, oo] is called a rate 
function if I is not identical oo and has compact level sets, i.e., if I~^([0, m]) = {x G A" : I(a;) < m} is compact 
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for any to > 0. A sequence {Xi}^J^ of A"- valued random variables Xi or the corresponding sequence {P;}^^ 
of probability measures on the Borel subsets of X satisfy the large deviations upper bound with speed a/ and 
rate function I if, for any closed subset C X , 

limsup — InP, [Xi G C) = limsup — InP, (C) < - inf I (a;) , (5.22) 
and they satisfy the large deviations lower bound if, for any open subset Q oi X , 

lim inf — In P/ {Xi £ G) = lim sup — In P/ (g) < - inf I (x) . (5.23) 

i^+co ai i^oo ai Q 

If both, upper and lower bound, are satisfied, one says that {Xi}^^ or {P;}^^ satisfy a large deviations 
principle. The principle is called weak if the upper bound in (|5.22p holds only for compact sets C. This notion 
easily extends to the situation where the distribution of Xi is not normalized, but a sub-probability distribution 
only. Observe also that one of the most important conclusions from a large deviations principle is Varadhan's 
Lemma, which says that, for any bounded and continuous function (/3 : A" ^ R, 

lim — In / exp {anp (Xi)) dP = - inf {I (x) - ip (x)} . 

i— ►+00 ai J x£X 

For a comprehensive treatment of the theory of large deviations, see |19) . 
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